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ABSTRACT 

Wc study the formation of low-mass and extremely metal-poor stars in the early universe. Our 
study is motivated by the recent discovery of a low- mass (M,, < -8 M ^,) and extremely metal-poor 
(Z < 4.5 X 10^^ Zq) star in the Galactic halo bv ICaffau et al.l (|2011f ). We propose a model that 
early supernova (SN) explosions trigger the formation of low-mass stars via shell fragmentation. We 
first perform one-dimensional hydrodynamic simulations of the evolution of an early SN remnant. 
We show that the shocked shell undergoes efficient radiative cooling and then becomes gravitationally 
unstable to fragment and collapse in about a million years. We then follow the thermal evolution of the 
collapsing fragments using a one-zone code. Our one-zone calculation treats chemistry and radiative 
cooling self-consistently in low-metallicity gas. The collapsing gas cloud evolves roughly isothermally, 
until it cools rapidly by dust continuum emission at the density 10^^-10^^ cm~^. The cloud core 
then becomes unstable and fragments again. We argue th at early SNe c an tri gger the formation of 
low-mass stars in the extremely metal-poor environment as ICaffau et al.l (|2011f ) discovered recently. 
Subject headings: stars: formation — supernovae: general 



1. INTRODUCTION 

Stars with masses M* < Q.SMq formed in the early 
universe may survive to the present day. Such old low- 
mass stars should typically have very low metallicites 
because they were formed from a gas that had not 
been significantly enriched with heavy elements. They 
carry invaluable inf ormation on the early chem ical evolu- 
tion of the Galaxy (jBeers fc Christliebll2005l) . Recently, 
ICaffau et al.l (I2011D discovered a low-mass and extremely 
metal-deficient star in the halo of the Galaxy. The star 
has extremely small abundances of carbon and nitro- 
gen. The upper limit of corresponding r netallicity of 
the s tar is estimated to be 4.5 x 10~^ Zq (jCaffau et al.l 
1201 It) . Furthermore, the ratio of the abundances of car- 
bon and oxygen derived from the spectra of the star 
suggests that the parent cloud from which the star was 
born had been chemically en riched by core-collapse su- 
pernovae (jCaffau et al.ll20li[) . 

It is generally thought that metal-free primordial stars 
are predominantly massive (e.g., Bromm et al. 2Qp9l) a l- 
though not extremely massive ([Hosokawa et al.l |2011[ ). 
This is largely because metal-free gas lacks efficient 
coolants other than hydrogen molecules; a primordial gas 
cloud is thermally stable throughout its pre-stellar col- 
lapse (Omukai & Yoshii 2003; Yoshida et al. 2008). The- 
oretical studies suggest that a trace amount of heavy el- 
ements provide efficient cooling in low-metallicity gas, 
whi ch may then enable the f ormation of low-mass stars 
(see iBromm fc Yoshidall2011l for a review) . There is an- 
other argument that the existence of dust grains in low- 
metallicity gas is cruci a l for the low-mass star form ation 
(iSchneider et al.l [200l lOmukai et al.l [20051 IMol) . Al- 
though it is still unclear whether or not gas metallicity 
is the key quantity to determine the characteristic stellar 
mass, and whether or not there is "critical metallicity" 
for low-mass star formation, there must be a mechanism 

^ E-mail: gen. chiakiOutap.phys.s.u-tokyo. ac.jp 



for the formation of low-mass stars in extremely low- 
metallicity gas^. 

SN explosions of the first generation stars are thought 
to have enriched the early universe with the first heavy 
elements. Three-dimensional sim ulations of ea r ly SN 
expl osions are perforin ed by, e.g., iBromm et al.l ()2003[ ) 
and iGreif et all ()2010[ ). These studies show heavy el- 
ements a re effectively dispers ed into the inter-galactic 
medium. IWhalen et al.l (|2010f ) argue that ultraviolet ra- 
diation from a massive progenitor star modifies the den- 
sity structure of the ambient gas before SN explosion. 
Overall, an early SN explosion can trigg er the formation 
of low -mass and low-metallicity stars. iNagakura et al.l 
()2009t ) study the evolution of SN remnants with metal- 
ficities of Z = 10"'' and IQ-^ Zq. The SN shell becomes 
gravitationally unstable for a certain range of parame- 
ters such as explosion energy, ambient gas density, and 
metallicity. They show that the typical mass scale of 
the fragments is ~ lO'^-lO^ Mq and argue that a low- 
mass s tar forms after the fu rther evolution of the frag- 
ment. iMachida et al.l (|2005| ) semi- analytically calculate 
the evolution of SN shells before and after the shell frag- 
mentation in a primordial gas. They show that the Jeans 
mass of the fragment is below the solar mass when the 
fragment becomes optically thick. However, it remains 
uncertain whether the gas cloud further fragment or not. 

In this paper, we study SN-triggered star forma- 
tion, using hydrodynamic simulations. We first run 
one-dimensional hydrodynamic simulations to follow the 
thermal evolution of an early SN remnant. We here as- 
sume the ambient gas around the SN has nearly the pri- 
mordial composition with a trace amount of heavy ele- 
ments, Z = 1.0 X 1O~^Z0 and 4.5 x lO'^Z©, presum- 
ably ejected by earlier generation stars. We determine 
whether or not the shocked SN shell is gravitationally 

^ In this paper, we use the term "low-mass" for a mass less than 
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Fig. 1. — We plot the evolution of the density, temperature, pressure, and velocity profiles of the supernova remnant for (EgN, no, Z) = 
(1 X 10^2 (,rg, 100 cm-3, 1.0 X 10"^ Zq) at the time (1) t = 3.2 X 10^, (2) 3.2 X 10"', (3) 3.2 X 10^, and (4) 8.6 X 10^^ yr. The blast wave 
driven by SN explosion sweeps up the ambient medium and the SN-shocked shell is formed. Epoch (1) corresponds to the Scdov- Taylor 
phase, in which the shell expands adiabatically. Then, the cold dense shell forms by radiative cooling at epoch (2), and eventually the shell 
is decelerated sufficiently to be gravitationally unstable at epoch (4). 



unstable to fragment by using linear analysis of density 
perturbations. We show that gravitational instability is 
triggered in the shell for a wide range of the supernova 
explosion energies, E^^ = 1 x 10^^-3 x 10^^ erg, and 
for large ambient gas density, ng > 10 cm^"^. Finally, 
we study the subsequent evolution of the collapsing frag- 
ments, using a one-zone code. While the gas cloud con- 
tinues collapsing gravitationally, it cools rapidly by dust 
thermal emission and eventually fragments to sub-solar 
mass clumps. 

2. NUMERICAL METHOD 

2.1. Code 

2.1.1. Ante- fragmentation phase 

We perform one-dimensional hydrodynamic calcula- 
tions to follow the dynamics of an SN remnant un- 
til the shocked gas shell undergoes gravitational insta- 
bility (hereafter, we call th is phase ante-fragment ation 
phase). We use the code oflKitavama et al] ()2004D and 
iKitavama fc Yoshidal (|2005[ ) which employs the second- 
order Lagrangian finite-difference scheme in spherically 
symmetric geometries. The basic equations are presented 
in §2 of Kitayama et al. (2004). 

We implement the relevant radiative cooling processes 
to the hydrodynamical code as follows. In the high tem- 
perature regime (T > 10"* K), assuming that the gas is 
in collisional ionization equ ilibrium, we ad opt the cool- 
ing function of Sutherland fc Dopital (I1993D . The cooling 
functions for Z < 10~^ Zq are derived by interpolating 
the cooling functions for Z = 10"'^ Zq and Z = 0. The 
cooling rate due to inverse Compton scattering is also 
included. 

In the lower temperature regime (T < lO** K), we 
solve the non-equilibrium chemical reactions and radia- 
tive processes of the following fifteen species: e~, H, 11+ , 
H", D, D+, D", He, He+, He", Ha, H+, HD, HD+, 



and HeH+. T h e mol ecular co oling rates are taken from 
IGalh fc Pallal (11991 for H2, iLipovka et aP ((20051) and 
from lFlower et all (|2000f ) for HD. We also calculate the 



fine-structure transition line cooling of heavy elements 
including C, O, Si, and Fe. These species are assumed to 
be in the form of C II, O I, Si II, and Fe II, respectively. 
The fraction of each species in the gas is derived by mul- 
tiplying the solar abundance by metallicity (Z/Zq). The 
solar abundances of C, O, Si, and Fe relative to hydrogen 
arc 2.45 x 10 "^, 4.90 x IQ-^, 3 . 63 x I Q-^, and 3.16 x 10-^ 
respectively ()Santoro fc Shull|[2006[ ) throughout the one- 
dimensional calculati on. Also, the atomic lev el transition 
rates are taken from lSantoro fc Shulll ()2006( ). 

2.1.2. Postfragnientation phase 

After the supernova shell becomes gravitationally un- 
stable (sec t j3. 21 below), fragments in a shell are expected 
to collapse further to trigger star formation. We perform 
one-zone calculations to follow the thermal and chemical 
evolution of one of such collapsing gas clouds. We assume 
that the collapse occurs approximately in a spherically 
symmetric manner. Then, the time evolution of the gas 
density, p, is described as 



dp 
'dt 



P_ 

ts' 



In this equation, tg denotes the free-fall time as 
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where G is the Newtonian gravitational constant. The 
evolution of specific internal energy, e, is 
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Fig. 2. — Wc plot the evolution of the density and tempera- 
ture of a fluid element, which is initially located at 20.8 pc from 
the center. Here, the result for the SN model of (_EsN,no,^) = 
(lxlO'^2 erg, 100 cm-3, 1.0x10-5 Z©) is shown. Solid and dashed 
curves denote the evolution before and after fragmentation, respec- 
tively. Both results are obtained by our one-dimensional calcula- 
tions. Points (1) to (4) indicate the epochs in Figuro[T] After epoch 
(4), the density in the region turns to increase by self-gravity. 

where p denotes the pressure and T and A denote heat- 
ing and coohng rate per unit volume, respectively. The 
temperature of the cloud, T, is related to the internal 
energy as 

e= 7 , (4) 

7-1 /imn 

where ^ is the mean molecular weight, ttih is the mass of 
hydrogen nuclei. In equation ((3]), 7 is the adiabatic index 
whic h is calculated by takin g average over all species as 
e.o.. lOmukai fc Nishil (ITqQSI) . 

In our one-zone model, the low-metallicity gas is as- 
sumed to consist of the following chemical species: e^, 
H, H+, H-, D, D+, D-, He, He+, He", C, C+, O, 0+, 
H2, R^, H+, HD, HD+, HD-, OH, H20,H20+, O2, Oj, 
CO, C0+, CO2, and COj. We assume that 72 % of car- 
bon and 46 % of oxygen are co ndensed to dust gra ins as 
in the local interstellar matter (jPoUack et al.|[l9M) . We 
include atomic cooling by C II and O I but neglect the 
contributions from less abundant Si II and Fe II because 
these species do not have effect on the thermal evolution 
of the gas cloud ([Schneider et al. 2012) ■ The adiabatic 
comprcssional heating and the chemical heating by H2 
formation are the dominant heating processes. We also 
consider the dust cooling in the calculations. The reac- 
tion rates and corresponding cooling rates are taken from 
[Onnikai et al, (2005,, 2Qi3)- We use the Sobole y length 
appro xima tion for optical l y thic k cooling as in lOmukail 
((2000h and lYoshidaeFaH mM) . 

2.2. Initial Conditions 

We set metallicities of the circumstellar medium Z — 
1.0 X 10^^ Zq conservatively and 4.5 x 1 0"^ Z(^, the 
upper bound of the metallicity of the star ICaffau et al.l 
(|2011l) observed. We assume that the progenitor star 
is hosted by an early small dark matter halo with mass 
Mh = 10^ Mq at a redshift of z = 20 (jYoshida et al.l 
|2003[ ). The gravitational potenti a l of th e dark matter is 
fixed as in iKitavama &: Yoshidal (|2005[ ) . We work with 
the ACDM model with the matter density fim = 0.27, 
the baryon density fib = 0.045, the cosmological constant 
J7a = 0.73, the Hubble constant h ~ 0.71, in units of 
100 km s"^ Mpc"^, and the absolute density fluctuation 



Of, = 0.80 consistent with WMAP 7 data ([Larson et al.l 

[mil). 

The density of the circumstellar medium at the time 
when SN explosions occur depends o n a number of factors 
such as the mass of the host ing halo ([Whalen et al.l[2003 . 
l2008t IKitavama et al.ll2004D and the ultraviolet luminos- 
ity of the progenitor star (IShu et al.|[2002[ : lAlvarez et al.l 
|2006() . Thus, for simplicity, we assume that the sur- 
rounding gas density is initially uniform for our one- 
dimensional hydrodynamic simulations and the initial 
density is set no =1, 10, 100, and 1000 cm^"^. 

We assume that the surrounding gas with primor- 
dial composition is pre-enriched with a trace amount of 
heavy elements. Thus, through this paper, the hydrogen 
mass fract ion is X = 0.74 , the helium mass fraction is 
Y = 0.26 ([Aver et al.ll2010[ ). an d the deuterium n umber 
fraction is D/H = 2.9 x 10"^ ([locco et al.ll2009t) . We 
set the initial gas temperature T = 10'^ K, and initial 
ionization fraction of hydrogen is uniformly 10"^. The 
fraction of hydrogen molecule is lO^'^ at the beginning^. 
The radial velocity is set v = km s^^. We test three 
explosion energies, i?sN = 1 x 10^^, 1 x 10^^, and 3 x 10^^ 
erg. We assume a point-source explosion and at the time 
t ~ 0^ the corresponding thermal energy is deposited in 
the innermost cell. 

3. RESULTS 

3.1. Shell Evolution 

3.1.1. Larger ambient gas density 

Figure [1] shows the structures of the shell for 
{EsN,no,Z) = (lxl0^2erg,100cm-3,1.0xlO-^Zo) at 
the time (1) t = 3.2 x 10^, (2) 3.2 x 10^, (3) 3.2 x 10^ and 
(4) 8.6 X 10^ yr. Figure [2| illustrates the evolution of the 
gas density and temperature for a fluid element that is 
densest when the shell becomes gravitationally unstable 
and that is initially located at 20.8 pc from the center. 
The evolution of the SN shell structure is described as 
follows. In the initial Sedov- Taylor phase, the SN rem- 
nant expands adiabatically as long as the radiative cool- 
ing time, tcooi = pe/A, is longer than t he shell expansion 
time scale, <cxp = R/V ([Sedovl[l958r ). Then, radiative 
cooling becomes effective and the gas just behind the 
shock starts to lose its thermal energy. At < ~ 2 x 10^ 
yr after the SN explosion, the cold shell is pushed by the 
large pressure of the inner hot gas, forming a very dense 
shell (the pressure-driven snowplough phase) . During the 
pressure-driven phase, the fluid element plotted in Figure 
[Renters the shock front, increasing its density and tem- 
perature rapidly. The inner hot region still impacts the 
cold shell, generating multiple weak shocks behind the 
shell. We sec this as oscillatory features from epoch (3) 
to epoch (4) in Figure [2| Finally, at epoch (4), the shell 
is decelerated sufficiently to be gravitationally unstable 
as we shall discuss in §3.21 

3.1.2. Smaller ambient gas density 



^ The structure and evolution of SN shells arc indeed insen- 
sitive to the circumstellar temperature and chemical compositions 
because the upstream pressure is negligible against the downstream 
pressure at the beginning of the shell evolution. Wc also confirm 
that at the end of simulations the upstream temperature becomes 
below 10^ K by radiative cooling even if we start the simulation 
with the temperature of 10* K. 
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Fig. 3. — Same as Figure 1, but for {Esn,riQ,Z) 
(3) 3.2 X 10*^, and (4) 1.0 X 10^ yr. 

Figure [3] shows the structures of the shell for 
(£;sN,no,^) = (1 X 10^1 erg, 1 cm-^^ i.o x lO^^ Zq) at 
the time (1) f = 1.0 x 10^ (2) 5.6 x 10^ (3) 3.2 x 10^, 
and (4) 1.0 x 10^ yr. Figure |4] illustrates the evolution of 
the gas density and temperature for a fluid element ini- 
tially located at 34.3 pc from the center. In this model, 
the shell is in the Sedov- Taylor phase around epoch (1) 
and then enters the pressure-driven snowplough phase 
at t ~ 2 X 10^ yr. In this case with low ambient gas 
density, the transition from the Sedov- Taylor phase to 
the pressure-driven phase occurs later than the model 
with high ambient gas density. When the shell sweeps 
up denser ambient medium, radiative cooling of the shell, 
which is generally proportional to the square of gas den- 
sity, becomes larger, and the enhanced radiative cool- 
ing makes the shell evolve faster. The shell continues to 
expand even after the pressure inside the shell declines 
due to radiative cooling, conserving the momentum (the 
momentum-conserving snowplough phase) . In our calcu- 
lation, the shell enters the momentum-conserving phase 
by epoch (3). The shell thickness increases because of 
the pressure in the shell. The internal energy of the 
shell remains higher by less effective radiative cooling. 
Around epoch (4) , the ellapsed time reaches the dynam- 
ical timcscalc for the dark matter halo hosting the pro- 
genitor star (idyn ^ 10^ yr for 10^ Mq halo at z ^ 20). 

3.2. Shell Fragmentation 

3.2.1. Criteria j or fragmentation 

There are several criteria proposed for the instability 
of an expanding and decelerating gas shell. lElmegreenI 
(|1994| ) derive approximate conditions for the collapse of 
an expanding and decelerating shell by linear analysis 
of density p erturbations under t he th in-sh ell approxima- 
tion. T hen, llwasaki &: Tsurib^ (|2008[ ) and llwasaki et al.l 
(|2011al lbl) include the eff ect of the internal structure of 
the shell to the result bv lElmegreenI (|1994[ ). For a shell 
with average column density over the all solid angle, Eq, 



(1 X 10^^ erg, 1 cm-^*, 1.0 X 10"^ Zq) at the time (1) t = 1.0 X 10^, (2) 5.6 X 10^, 



shell radius, i?, and velocity y, the angular wavenumber 
of the perturbations that grow most rapidly is 



^pcak 



~5 — ^' 



(5) 



and the growth rate of the perturbations with wavenum- 
ber r^poak is 



ijj- 
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In equations ([S]) and © , Coh is the effective sound speed 
as 



Ceff 
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where p,nax is the peak density in the shell and ct = 
{kT/iiniiiy^^ is the sound speed for an isothermal shell 
with temperature T. The numerical factor A is deter- 
mined as 0.39 by the detailed linear analysis for a shell 
with its outer boundary defined by the shock front and 
with its inne r boundary defined by t he stellar wind from a 
central star (jlwasaki et al.ll2011al lb[) . Note that replacing 
Coff with Ct and Eg with pqR/3, where po is the average 
pre-shell ma ss density, i i i equa tions ([5]) and (O recovers 
the result of lElmegreenI ()1994D . 

In equation (|6]), the first and second terms in the 
right hand side represent the effect of the shell expan- 
sion and the last term represents the effect of the shell 
internal pressure and self-gravity. If the peak growth 
rate, Wpeak, is positive, perturbations with the corre- 
sponding wavenumber have a growing solution, which 
means the SN-driven shell becomes gravitationally un- 
stable and fragments with scales corresponding to J^pcak 
start to grow in the shell. It is obvious that fast expan- 
sion (large V) tends to smooth out perturbations. For 
the gas to become gravitationally unstable, its expansion 
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Fig. 4. — Same as Figure 2, but for {Esj<i,no,Z) = (1 x 
10^^ erg, 1 cm~^,1.0 X lO"'' Zq) . Points (1) to (4) indicate 
the epochs in Figure [3] Dotted line shows an adiabatic trajec- 
tory, T on nri~^. The trajectory is along an adiabatic line between 
epochs (2) and (4) because the shell expands radially for small ex- 
ternal pressure, while in larger ambient density case, the density 
declines slowly because the shell expansion is suppressed by the 
large external pressure and the temperature decreases by radiative 
cooling (Figure [2J. 

velocity must be sufficiently small such that 



V ■kGY.q 



3.2.2. Application of the criterion to our results 



(8) 



In this section, we discuss whether the gravitational 
ins tability can be t riggere d or not. We apply the criterion 
by llwasaki et al.l ()20f fal fo) to our calculations because 
the internal structure of a shell is taken into account in 
their analysis^. In order to investigate the stability of the 
shell, we calculate the time evolution of the peak growth 
rate, Wpoak [equation ([6])] by using the radius of the shock 
front, the shock velocity, volume- weighted mean density, 
mass-weighted mean temperature, and mean molecular 
weight in the shell. In the snowplough phase, the shell 
is defined by the outer boundary (shock front) and the 
inner boundary. The outer boundary is determined as 
position where the density is 1.1 po- The displacement of 
the inner boundary is defined as the point where density 
recovers the pre-shell density, poj inside the shock radius. 

After the shell becomes gravitationally unstable, the 
shell is assumed to be separated into fragments. If the 
time for fragmentation is comparable to or larger than 
the dynamical timescale of the hosting halo, idyn, we do 
not count such cases as successful fragmentation of the 
SN shell in the cosmological context. After the dynami- 
cal timescale, it is likely that the halo itself collapses or 
is incorporated into a large halo, thereby the structure 
of the entire system would be much different. 

Applying the fragmentation condition, we find 
that for the SN model of {Eo,^,nQ,Z) = (1 x 
10^^ erg, 100 cm~^, 1.0 x 10~^ Z©), the compressed and 
decelerated shell becomes gravitationally unstable at 



* llwasaki et al.l l|2011al lH) presented the linear analysis for a 
shell defined by a shock front and an inner boundary formed by 
the pressure of the inner region. It is thus more appropriate to 
our results of the SN-driven shell defined by the shock front and 
th e inner boundary alth ough the shell inner boundary in the study 
of llwasaki et al.1 I l2011al lbl) is driven by the pressure of the stellar 
wind from the central star. We have also checke d that the r e sult i s 
qualitatively unchanged if we apply the result bv lElmegreenI 1)19941 ') 
which ignores the effect of the shell structure. 



iinst — 8.3 X 10^ yr. By contrast, for {Ea,^,nQ,Z) = 
(1 X 10^^ erg, 1 cm~^, 1.0 x 10~^ Z©), the peak growth 
rate Wpoak becomes positive at tinst = 1.0 x 10^ yr when 
the shell is about to fall back {V ~ 0) by the gravita- 
tional potential of the host halo. In this model, since the 
time iinst is comparable to the dynamical timescale, it is 
uncertain whether the shell fragments or collapses into 
the center of the halo. 

Table[T]shows the physical properties of the shell at the 
time tinst for SN models we study. The results are shown 
only for Z = 1.0 x 10^^ Zq because they are almost ex- 
actly the same as the results of the Z = 4.5 x 10"^ Zq 
models. It is because the fine-structure cooling of met- 
als such as C and O has little effect on the evolution 
of the SN remnant in such an extremely metal-poor 
regime until the shell becomes gravitationally unstable. 
Also, the time iinst is insensitive to the SN explosion 
energy. For a radiative remnant, the shock radius in- 

1/7 

creases only as i? oc EA -j in the pressure-driven phase 
(jOstriker fc McKeel[T988l) . The time tjnst depends on the 
ambient gas density. As the ambient gas density is larger, 
denser and colder shell is formed as we discuss in ij3.1.2l 
and the growth rate of the perturbations on the shell, 
i^pcak becomes positive earlier. The Jeans mass (typi- 
cal mass) of the fragment is shown in the last column 
of Table [H The mass decreases as ambient gas density 
increases because the density of the shell swept up by the 
blast wave becomes larger. 

Figure [5] shows the results of all the models that 
we investigated in a (E'snj'^o) plane. These results 
are essentially the same for the respective model with 
4.5 X 10^^ Zq. The circles indicate models where both 
our criteria are met. For models with ambient density 
no = 1 cm~'^, the time for the gravitational instability, 
<inst, is comparable to tdyn- These models are marked by 
triangles. We cannot determine from our ID calculations 
whether the shell fragments before the host halo evolves 
significantly. It would be needed to run a cosmologi- 
cal simulation that incorporate the hierarchical growth 
of dark halos, in order to follow the evolution of such 
systems properly. 

3.3. Evolution of Fragments 

We investigate the subsequent evolution of the frag- 
ment for the models with large ambient gas density 
no > 10 cm~"^ where the SN-driven shell becomes grav- 
itationally unstable well before 10^ yr. In these models, 
the time, iinst when SN shells fragment is comparable to 
the free-fall timescale, tg, of the shell at the time tinst- 
The fragments are expected to contract gravitationally 
on the free-fall time scale (equation [ij . Curves in Fig- 
ure ini denote the evolutionary trajectories of the core of 
each fragment calculated by our one-zone code for the 
metaUicities 1.0 x 10"^ Z© and 4.5 x 10"^ Zq when SN 
explosion energy and ambient density are (i?sN,'T-o) = 
(1 X 10^^ erg, 100 cm^'^). The initial density, temperature 
and the chemical species' abundances are taken from the 
result of our one-dimensional calculation. We use the 
result at the time when the gravitational instability is 
triggered in the shell. 

The temperature stays between a few hundred to one 
thousand Kelvins while the density increases over many 
orders of magnitudes. The compressional heating is 
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TABLE 1 
Properties of Fragments (Z : 



1.0 X 10" 



Zq) 



no 
(cm-3) 


-EsN 
(1051 erg) 


*inat 

(Myr) 


riH 
(cm-3) 


T 
(K) 


R 

(pc) 


(pc) 


^pcak 


{Mq) 


A^Joans 

(Mo) 


1 


1.0 
10.0 
30.0 


10.23 
10.84 
10.35 


2.9E+00 
6.7E+00 
l.lE+01 


568.6 
466.8 
419.7 


115.9 
196.6 
250.9 


22.84 

12.89 

9.17 


3.0 

8.7 
14.4 


3.0E+05 
1.3E+06 
2.6E+06 


4.8E+05 
2.4E+05 
1.6E+05 


10 


1.0 
10.0 
30.0 


3.39 
2.99 
2.63 


5.9E+01 
1.7E+02 
2.4E+02 


300.2 
239.8 
215.2 


45.6 
75.2 
93.9 


3.53 
1.72 
1.42 


6.9 
21.5 
36.5 


1.7E+05 
6.7E+05 
1.3E+06 


4.1E+04 
1.7E+04 
1.2E+04 


100 


1.0 
10.0 
30.0 


1.06 
0.83 
0.72 


l.lE+03 
2.4E+03 
3.8E+03 


191.6 
182.8 
186.4 


17.9 
28.5 
35.2 


0.68 
0.43 
0.33 


14.2 
38.5 

58.4 


9.3E+04 
3.5E+05 
6.4E+05 


4.9EH-03 
3.0E+03 
2.5E+03 


1000 


1.0 
10.0 
30.0 


0.32 
0.25 
0.22 


1.4E+04 
3.1E+04 
4.0E+04 


202.5 
206.5 
223.0 


6.9 
10.7 
13.2 


0.19 
0.12 
0.12 


19.5 
48.8 
69.8 


5.1E+04 
1.8E+05 
3.4E+05 


1.5E+03 
l.OE+03 
l.OE+03 



Note. — Z: mctallicity of the ambient material, no: ambient gas density, Eq^: explosion 
energy, finst^ time when the gravitational instability first appears, njj: volume-averaged Hydro- 
gen number density in the shell, T: mass- weighted average of temperature in the shell, R: radius 
of the shock front. A/?: shell thickness, 77poak: peak wavenumber of the perturbations, Msholi: 
mass of the shell, and Af jeans : Jeans mass, at t = tjnst 



largely balanced by radiative cooling until the Hydro- 
gen atomic number density increases to nn ^ lO^^cm"^. 
There are, however, two slight declines as clearly seen in 
Figure [S] The first decline appears at nn ~ 10^ cm~'^ 
where cooling by HD molecular lines is efficient. The sec- 
ond decline appears at nn ~ 10^^ cm~'^ where cooling by 
dust thermal emission is efficient. In the latter regime, 
the gas cloud core becomes further unstable a nd is likely 
to fragment into multiple clumps (jOmukai et a l. 2005). 
The gravitational instabilities occur when the dust 
cooling timescale is shorter than the free-fall timescale. 
The gas cooling by dust is due to heat transfer from gas 
particles to cold dust grains (see IStahler fc Fallal I2005L 
§7.4.3). Then, essentially the gas pressure determines 
the juncture for efficient dust cooling as 



p = nukT > 



32G fadPd 



9n V Z^a 



(9) 



where a^ and pd denote the mean radius and mass den- 
sity of an individual dust grain, Zd is dust mass frac- 
tion, and a is the accommodation coefficient for heat 
transfer between gas and dust. Fro m the typical va l- 
ues of ad ~ 10~^ cm, p^ — 3 g cm~^ (jLove et al.lll994[ ). 
Zd ^ 10~'^{Z/Zq), and a ^ 0.1, the region on the n-T 
plane where dust cooling becomes efficient is 



P 



nuT > 4.9 X 10^^ cm 



V 10-5 



(10) 



The triangles in Figurc[6]dcnotes the points where the gas 
clouds satisfy the condition given by equation ^TU\\ . The 
corresponding Jeans mass at the point is 0.4 Mq for Z = 
1.0 X 10-5 Zq and 3 Mq for 4.5 x lO'^ Zq. The Jeans 
mass further decreases as the gas temperature rapidly 
decreases by dust cooling, as can be seen in Figure [SI 

The rapid dust cooling causes defor mation of the cloud 
core into a filamentary structure. iTsuribe fc Omukail 
()2Q08[ ) show, using three-dimensional simulations, that 
an initially perturbed gas cloud is significantly elon- 
gated to filamentary shape and further fragments into 
small pieces when dust cooling operates for the nictal- 
hcities 10"^ < Z/Zq < lO'^. The deformation in- 
stability occurs when the deformation parameter, de- 
fined as the effective ratio of specific heats, 7cff = 



dlogjp/dlog/o, is less than a critical va lue 7crit — 1-097 
(jHanawa fc Matsumotoll2000l: lLaj||2000[ ). In our calcula- 
tions, 7eff declines rapidly below the critical value when 
the dust cooling becomes efficient. The deformation pa- 
rameter is 0.5 at a minimum when nn ~ 10^^ cm"'^ for 
Z = 1.0 x 10-5 Zq. The value is common with all ex- 
plosion energies we set. For Z = 4.5 x 10^5 Zq, the 
minimum value of 7off is 0.3-0.4 at njj ^ 10^^ cm-^. 
The contracting gas clouds are strongly unstable to de- 
formation in these regimes. 

When the instability occurs, regardless of the SN ex- 
plosion energies and metallicities, the Jeans mass is 
0.01-0.1 AIq which is the characteristic mass of a frag- 
ment (protostar). The ambient gas quickly accretes 
onto the proto stars and their mass likely become larger 
(|Larson|[l969D . However, the subsequent evolution of 
the protostars will be highl y dynamical. According to 
iPortegies Zwart et al.l ()2010l) , some of the protostars can 
attain velocities above the escape velocity of the parent 
gas cloud through dynamical interactions among them. 
They would then escape from the gas cloud before the 
protostars become massive. One natural outcome would 
then be the formation of low-mass stars. 

4. DISCUSSION 

We have studied the evolution of an early SN remnant 
with a mctallicity of ^ 10-5 Zq. The supernova shell sat- 
isfies the criterion for gravitational instability at i ~ 105- 
10^ yr, to form multiple fragments with the characteristic 
Jeans mass of Mjcans ^ lO'^-lO'* M© for a wide range of 
the SN explosion energies, Egjsi = 1 x 105^-3 x 105^ erg, 
and large ambient gas density, tiq > 10 cm-'^. We have 
then followed the thermal evolution of one of the frag- 
ments using one-zone calculations. The gas cloud core 
evolves roughly isothermally until the density increases 
to riH ~ 10^^ cm-'^, where efficient cooling by dust ther- 
mal emission brings the gas temperature from ^ 10"^ K 
to ~ 300 K rapidly. The unstable core is expected to 
deform and fragment again into multiple clumps (pro- 
tostars) with masses of ~ 0.01-0.1 Mq regardless the 
SN explosion energies and metallicities. If some of the 
protostars escape from the host gas cloud before grow- 
ing in mass over ~ 0.8 Mq, they would live long to the 
present day. We argue that the shell fragmentation of an 
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Fig. 5. — Fragmentation conditions for each our SN model. The 
result is common between the model with Z = 1.0 X 10~^ Zq and 
4.5 X 10~^ Zq. The circles indicate models where fragmentation 
condition is met. The triangles indicate models for which it is 
uncertain by our calculations whether the shell fragments or not 
(see text). 

early SN remnant is a viable scenario for the formation 
of low-mass and extremely low-metallicity stars. 

There have been several theoretical studies which sug- 
gest the formation of low-mass stars in low-metallicity 
gas using numerical simulations. For metallicites of 
r^ 10~^Zq, sma ll mass clumps are forrn ed from a ro- 
tating ga s cloud (iTsuribe fc Omukaill2008f ) and in turbu- 
lent gas (jDopcke et all 1201 1[ ) in manv conditions these 
authors examined. Our hydrodynamic calculations are 
based on SN explosions of an early generation of stars. 
The end point of the hydrodynamic simulation of the SN 
remnant, when the shocked shell becomes gravitationally 
unstable, provides physically motivated initial conditions 
for our one-zone calculations. Therefore, we are able to 
propose, for the first time, a consistent model of how 
low-metallicity stars are formed in the early universe. 

We note that there are still limitations in the one- 
dimensional and one-zone simulations. Non-spherical 
structures of an SN remnant mig ht have substantial 
effects on the shell fragmentation. | Chevalier fc Thevsl 
(|1975f ) and lChevalier fc Imamural ()1982t ) show that shell 
instabilities might enhance density perturbations around 
an SN remnant to form dense clumps in early epochs of 
the remnant evolution. Dense gas clumps could form 
even earlier than we have shown in this paper. Clearly, 
full three-dimensional calculations arc needed to study 
in detail the evolution of an SN shell. Our calculations 
may provide basic configurations for such studies on the 
formation of low-mass stars in low-metallicity gas. 

Currently standard models of galaxy formation are 
based on the hierarchical growth of cold dark matter 
halos. Recent high-resolution iV-body simulations show 
that tens of thousan ds of early sma l l halos assemble a 
Milky Way size halo (iGao et al.ll20Tnl:IWang et al.l[2nTlh . 
Although stars formed in the earliest halos are likely pop- 
ulated in the central bulge region, there should be a sub- 
stantial fraction of the remnants of early halos orbiting 



in th e halo of the Galaxy (iTumlinsonI 120101 : iGao et akl 
|2010( ). ICaffau et al.l (|2011| ) may have found one of such 
"fossils". Future surveys of metal-poor stars will pro- 
vide the relative number fraction of extremely metal-poor 
stars, which will then constrain the overall formation effi- 
ciency of low-mass stars in the scenario we propose here. 
The authors thank K. Omukai, M. N. Machida, K. 
Iwasaki, T. Tsuribe, and colleagues at IPMU for fruit- 
ful discussions and helpful advice. We are grateful to 
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entists (S: 20674003, B: 21740139) by the Japan Society 
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10" 



^ 3 



Sio^ 



10' 



. 1 1 1 1 1 


/ 1 .■■ 1 .■ 1 


■ 1 ■/ : 


Mjeans=10Vo 


IO^Mq IO^Mq 10^ Mq 


10°Me ; 




..Z=1 .OXl..0'^ Zq / 


/ 10-^3^/- 




■■' /^..t=4.5x)&^ ZqV 


■4^ 


- ,.-\ /" 


/' / / / 


/j'O'^Mg' 


1 ■■■■', 1 ■■■■', 


1 ■■■', 1 ■/; 1 


1 



10-= 



lO'* 



10° 



10° 



10 



10 



10 



12 



14 



10 



16 



Hydrogen number density [cm" ] 



Fig. 6. — We plot the thermal evolution of the core of a collaps- 
ing fragment for (-BsNi'^o) = (1 x 10^^ erg, 100 cm~^). Red and 
orange curves denote the trajectories for Z = 1.0 X lO"'' Zq and 
4.5 X 10~^ Zq, respectively. The initial condition of the one-zone 
calculation is determined using the result of our one-dimensional 
calculation. The initial density and temperature are set as the av- 
eraged values over the shell at the time when the shell fragments. 
Triangles are plotted at the poin ts where the rapid dust cooling 
becomes dominant [equation mOI I is satisfied] for each metallicity. 
Thick solid line indicates the opacity limit characterized by rj = 1, 
where the cloud core is optically thick to continuum radiation (see 
e.g. Onmkai 2000). The instability that induces deformation of the 
fragment occurs when the temperature decreases by rapid cooling 
(7cff ^ !)• The fragment continues to deform, inducing further 
fragmentation until the temperature stops to decrease. The corre- 
sponding Jeans mass at the point is 0.01—0.1 Mq which is regarded 
as the characteristic mass of a protostar. 
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